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We examine the constraints on the luminosity-dependent density evolution model for the evolution 
of blazars given the observed spectrum of the diffuse gamma-ray background (DGRB), blazar source- 
count distribution, and the blazar spectral energy distribution sequence model, which relates the 
observed the blazar spectrum to its luminosity. We show that the DGRB observed by the Large 
Area Telescope (LAT) aboard the Fermi Gamma Ray Space Telescope can be produced entirely by 
gamma-ray emission from blazars and nonblazar active galactic nuclei, and that our blazar evolution 
model is consistent with and constrained by the spectrum of the DGRB and flux source-count 
distribution function of blazars observed by Fermi-LAT. Our results are consistent with previous 
work that used EGRET spectral data to forecast the Fermi-LAT DGRB. The model includes only 
three free parameters, and forecasts that > 95% of the flux from blazars will be resolved into point 
sources by Fermi-LAT with 5 years of observation, with a corresponding reduction of the flux in the 
DGRB by a factor of ~2 to 3 (95% confidence level), which has implications for the Fermi-LAT's 
sensitivity to dark matter annihilation photons. 
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I. INTRODUCTION 

The source of the extragalactic isotropic diffuse 
gamma-ray background (DGRB) has been an unsolved 
question in astrophysics for some time. In this paper, 
we show how the DGRB spectrum can be produced by a 
combination of blazar and nonblazar active galactic nu- 
clei (AGN) gamma-ray sources. We also show that the 
blazar flux source-count distribution function (dN/dF) 
is consistent with the full DGRB originating from these 
sources. Furthermore, we show how less-detailed models 
of the blazar contribution failed to be consistent with the 
DGRB. We explore how the implications for dark matter 
detection or constraints from the DGRB will evolve as 
the blazar sources of the DGRB are resolved. 

The DGRB was first discovered by the SAS 2 exper- 
iment in 1975, for gamma-ray emission in the range of 
35 to 300 MeV jl,. This background was seen at ener- 
gies up to 20 GeV by the EGRET Collaboration, and 
it was confirmed at these energies in the first-year data 
from the Large Area Telescope (LAT) aboard the Fermi 
Gamma-Ray Space Telescope [2H4|. The assumed extra- 
galactic source of the DGRB is determined by measuring 
the complete diffuse (unresolved) flux and then subtract- 
ing off a model to account for the background coming 
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from our Galaxy. This yields a measure of the flux com- 
ing from unresolved diffuse sources, presuming there is 
no minimal isotropic component from the Galaxy, e.g. 
dark matter annihilation or decay. The DGRB has been 
used to constrain dark matter annihilation in Galactic 
and extragalactic sources [IHl]. 

The most recent measurement of the DGRB was per- 
formed by the Fermi-LAT. In the Fermi-LAT Collabora- 
tion analysis, the gamma-ray intensity was measured in 
the range of 100 MeV to 100 GeV above 10° in Galactic 
latitude (|&| > 10°). The total flux is modeled by stacking 
the spectra of known sources with the cosmic-ray back- 
ground, the Galactic diffuse background, and the DGRB. 
This analysis gives a DGRB intensity that is roughly 25% 
of the total observed flux. The DGRB seen by the Fermi- 
LAT is consistent with a power law in energy with index 
2.41. This value for the DGRB is notably softer at high 
energies than was previously seen in the EGRET Collab- 
oration, which is partly due to an updated model of the 
diffuse Galactic emission in Ref. [3] (hereafter FSIO). 

A detailed spectral energy distribution (SED) sequence 
model of blazars can reproduce the DGRB [8j[9]. We ex- 
plore this model in this work. Many models have been 
proposed to explain the DGRB. It has been shown that 
emission from AGN can account for the diffuse back- 
ground from 10 keV to 100 MeV, but above that en- 
ergy, this model cannot account for the large gamma-ray 
flux [101 . Radiation from star- forming galaxies could ac- 
count for much of the DGRB up to 10 GeV, but this also 
cannot explain the high intensities observed at higher en- 
ergies [11| . Emission from millisecond pulsars has been 
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proposed as a source as well [H]. However, millisecond 
pulsars as a dominant source of the DGRB may be in- 
consistent with the lack of anisotropy in the DGRB p^ . 

Dark matter annihilation, both as a component of 
the extragalactic diffuse emission and as an unaccounted 
foreground from the Milky Way can contribute to the 
DGRB, but the fluxes from dark matter are expected to 
be lower than the DGRB flux and have a different spec- 
tral shape [3 [6]. However, measurements of the DGRB 
are one of the strongest ways to constrain dark matter 
annihilation [7]. If dark matter is a significant contrib- 
utor, it may be disentangled from astrophysical sources 
due to its angular correlation on the sky [M]. Pioneer- 
ing work proposed that blazars could account for all of 
the DGRB seen by the EGRET Collaboration [H]. The 
blazar class of AGN has been studied in depth as the 
origin of the DGRB at high energies [T5H26] . 

In Ref. [5] it was shown that the DGRB can be com- 
posed of blazars and nonblazar AGN in the luminosity- 
dependent density evolution (LDDE) SED blazar model. 
This model contains only three free parameters describ- 
ing the gamma-ray luminosity function (GLF) of blazars. 
We show that this model is consistent with producing 
the full DGRB spectrum as well as the blazar source- 
count distribution, dN/dF, of blazars as measured by 
Fermi-LAT. In addition, we constrain this model by these 
measurements and find parameters for which the model 
successfully reproduces these measurements. Note that 
both the source-count distribution dN/dF and DGRB 
spectrum are predicted by the model, and not an input 
to the model. 

Recent work by the Fermi-LAT Collaboration found 
that the DGRB could not be composed entirely by 
blazars [4] (hereafter FBIO). However, that work adopted 
an over-simplification of the blazar SED to be a single 
power-law (PL), independent of blazar luminosity, which 
is inconsistent with the observed spectral luminosity de- 
pendence seen in the SED sequence [27 . In contrast, in a 
separate paper, the Fermi-LAT Collaboration emphasizes 
the need for including departures from pure-PL behav- 
ior in blazar spectra when calculating the contribution 
of unresolved low-luminosity blazars to the DGRB [5S]. 
Incorporating the SED departure and its dependence on 
blazar luminosity evolution when modeling the DGRB is 
exactly the intent of the work presented here. 

Furthermore, the blazar model in FBIO lacks a physical 
evolution model for blazars. Instead of the source-count 
distribution resulting from the cosmological evolution of 
blazars, the source-count distribution is an input to the 
model, as a broken power-law with four free parameters. 
Note that even though the model in FBIO is simplistic, it 
contains more free parameters than the LDDE plus SED- 
sequence model explored here. In our approach there are 
three parameters in the adopted blazar model which de- 
scribe the relation between the GLF and x-ray luminos- 
ity function (XLF). Because the FBIO model employs a 
pure-PL luminosity-independent SED with a broken-PL 
source-count distribution, the conclusions of that work 



do not apply to the model examined here. Other pa- 
rameters in our work {e.g., the SED sequence and the 
low-energy nonblazar AGN model) are constrained by 
other observations and remain fixed in our blazar model 
analysis. Namely, the observational constraints on the 
SED sequence come from spectral population models of 
blazars as in j27l, and the nonblazar AGN spectrum is 
constrained by the hard x-ray luminosity function derived 
from HEAOl, ASCA, and Chandra x-ray AGN surveys 

[iniEg]. 

A recent paper by Malyshev and Hogg [SOj using the 
one-point probability distribution function (PDF) of the 
DGRB also concludes that blazars cannot constitute the 
total DGRB flux as measured by Fermi-LAT, when mod- 
eled as a pure-PL SED with a fixed dN/dF. However, 
this conclusion also only applies to the model which 
they consider, which adopt blazars as having pure-PL 
luminosity-independent SEDs, and not to the LDDE 
SED-sequence model examined here. 

Because observed blazars make up about 15% of the to- 
tal gamma-ray flux, unresolved blazars are a likely candi- 
date to make up the DGRB Blazars were the most 
numerous point-source objects observed by the EGRET 
Collaboration 'WT . Additionally, observed blazar spec- 
tra tend to follow a similar power law in energy as the 
DGRB. However, it is known that blazars have a lumi- 
nosity dependence to their spectral shape, which is in- 
corporated in the SED-sequence model [57], but ignored 
in the analysis of FBIO. 

Blazars are the combination of two classes of AGN: 
fiat-spectrum radio quasars (FSRQs) and BL Lacertae 
objects (BL Lacs). FSRQs are AGN that have spectral 
index a,. < 0.5 in the radio band and have radio emission 
lines with equivalent width greater than 5 A. BL Lacs 
have no strong absorption or emission features, and have 
equivalent widths less than 5 A j32j. Broadly speaking, 
blazars tend to have their bolometric luminosities domi- 
nated by the gamma-ray luminosity and have great vari- 
ability in that luminosity. Therefore, it is believed that 
blazars represent the small set of AGN that are observed 
along the jet axis, as opposed to nonblazar AGN which 
are observed far from the jet axis and dominate emis- 
sion by their luminous accretion disk. This jet source 
is expected to be relativistically beamed, as opposed to 
the more isotropic flux coming from the AGN's accretion 
disk [51 [33. 

Different models of blazar emission have been pro- 
posed in the literature [T5Vl26j . One is the pure lu- 
minosity evolution (PLE) model of the distribution of 
blazars HH UHl [Ml ■ In this model, only the blazar 
luminosity is evolved in redshift. An alternative model, 
LDDE, relates the gamma-ray luminosity of blazars to 
the redshift-dependent distribution of x-ray emission 
from nonblazar AGN [53]. This technique more realis- 
tically fits the blazar evolution to the AGN distribution, 
rather than assuming that all blazars have identical evo- 
lution regardless of luminosity. In many models for blazar 
spectra, a simple power-law or distribution of power laws 
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is used as the intrinsic blazar spectrum, but more detailed 
frequency-dependent models have been used as well |22j . 

Here, we employ the LDDE model for blazar distri- 
butions. For the intrinsic spectrum of blazars, we use a 
frequency-dependent SED based on the multiwavelength 
study of Ref. [21] ■ We use these models to derive the 
differential blazar spectrum in redshift, luminosity, and 
energy. By integrating over these variables, we can deter- 
mine the number of detectable blazars for given detector 
sensitivities, and we can calculate the expected gamma- 
ray flux from unobserved blazars to determine how sig- 
nificantly they contribute to the DGRB. Additionally, we 
add a fixed nonblazar AGN component to our predicted 
blazar flux, which should make the net flux from our 
model fit the diffuse background over the energy range 
from 10 keV to 100 GeV. 

Below, we begin by describing the DGRB seen by the 
Fermi-LAT as well as its data on blazars. We will then 
describe our model in detail, specifying the evolution 
model and SED used in our calculations and how we fit 
these to the known data. We use this model to predict 
the ability of the Fermi-LAT to detect blazars and how 
this will affect the DGRB. Throughout the paper, we take 
a flat universe with the cosmological parameters f2m = 
0.272, riA = 0.728, and Hq = 70.2 km s'^ M-pc'^ [§5]. 
Note, the use of h in the text refers to Planck's constant, 
and not the Hubble parameter. 



II. FIRST- YEAR FINDINGS BY THE 
FERMI-LAT COLLABORATION 

A. DGRB Measurements 

From its first year of data, the Fermi-LAT has mea- 
sured a spectrum for the DGRB (FSIO). To get this spec- 
trum, the total gamma-ray intensity had known sources 
subtracted from it, as well as the background from cosmic 
rays, and the expected Galactic diffuse emission. At this 
time, resolved extragalactic sources account for about 
15% of the total gamma-ray flux in the sky. To calculate 
the gamma-ray emission from Galactic cosmic rays, the 
local cosmic-ray spectra are extrapolated to give source 
populations, which are then propagated through appro- 
priate target distributions using the GALPROP particle 
propagation package [351 [S3 • This diffuse Galactic emis- 
sion is the largest component of the DGRB, comprising 
roughly half of the total observed intensity. A small com- 
ponent to the DGRB is a background due to cosmic-ray 
interaction with the Fermi-LAT itself. This background 
has been studied in detail in FSIO and is very well char- 
acterized. This background accounts for 1 to 10% of the 
total emission, with a greater fraction at low energies and 
a lesser fraction at high energies. The residual intensity 
after all of these components have been removed is called 
the isotropic DGRB. It makes up around 25% of the total 
emission. Because of the model dependence of these sub- 
tractions, the uncertainties on the DGRB are dominated 



by systematics (FSIO). The DGRB may come from un- 
resolved extragalactic sources or unaccounted Galactic 
sources, such as millisecond pulsars, or, potentially, from 
Galactic dark matter annihilation or decay. 



B. Point-Source Sensitivity 

The Fermi-LAT detector has a spectrally dependent 
point-source sensitivity due to the higher spatial resolu- 
tion of the instrument to higher-energy photons. The flux 
limit to point sources is shown in Fig. [l] along with the 
sample of blazar fluxes and spectral indices from FBIO. 
In FSIO, the DGRB spectrum is compared to that mea- 
sured by EGRET, which had a point-source sensitivity 
of 1 X 10~^ ph cm^^s"""^, despite the fact that the point- 
source sensitivity of the two instruments, and therefore 
the measured DGRB flux between the two instruments' 
measurements, are quantitatively different]^ We derive 
the flux limit from the sample of blazars used in FBIO, 
using the lowest-flux end of the blazar sample, which 
satisfied the test-statistic TS — 25. In FBIO, the source- 
count distribution and DGRB spectrum was fit with only 
blazars resolved at TS = 50; therefore, the point-source 
limit is augmented by a factor of 2, as shown by the 
solid in Fig. [T] with the point-source sensitivity always 
below or equal to Fermi-LAT's believed completeness for 
all spectra sources at 7 x 10~* ph cm~^s~^. 

Importantly, it should be made clear that a fixed 
point-source sensitivity cannot be exactly specified for 
the DGRB spectrum derived in FSIO. In that work, all 
sources above a TS = 200 are allowed to vary in the 
amplitude of their flux during the fitting of the extra- 
galactic isotropic DGRB. Therefore, the exact fiux-limit 
of the DGRB spectrum, and therefore the nature of the 
spectrum itself, as presented in FSIO, is ill-defined. We 
therefore adopt the best-estimate method of modeling 
the DGRB spectrum as done by the Fermi-LAT Col- 
laboration itself in FBIO, with a TS = 50 spectrally- 
dependent fiux limit. We define the power-law photon 
index F for the non-power-law SED-sequence model of a 
blazar by fitting a power law to the Poisson-limited spec- 
trum within the observed energy range of Fermi-LAT. 

As the point-source sensitivity of Fermi-LAT improves 
with integration time, the resolution of the extragalactic 
DGRB into point sources will not proceed proportionally 
to the sensitivity, but rather in a combination of the sen- 
sitivity with where the population of extragalactic emit- 
ters lies with respect to that sensitivity /spectral-index 



Because of this direct comparison in FSIO, in the vl preprint 
of this work, a point-source sensitivity cutoff of the mea- 
sured DGRB spectrum of FSIO was adopted to be 1 X 
10~^ ph cm~^s~^, instead of the spectrally dependent sensitiv- 
ity here. This does not change our conclusions, but does modify 
our best-fit model parameters and our 5-year forecast DGRB 
spectra. 
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FIG. 1. Shown is a sample of the blazar gamma-ray fluxes 
above 100 MeV (-Fioo) versus their power-law fit spectral- 
index r from FBIO. The blazars (points) are shown above 
point-source detection test-statistic TS — 25 (with the corre- 
sponding point-source limit shown as the dashed line), while 
those below TS = 50 are modeled, in our work and in FBIO, 
to contribute to the DGRB as measured by FSIO (with point- 
source limit shown by the solid line) . Note the total luminos- 
ity vs spectrum dependence of the blazar population evident 
in this plot. The first-year Fermi-LAT point-source sensitiv- 
ity is complete above the dashed line at 7 x 10~* ph cm~^s~^ 
(FBIO). 

plane. In particular, for the LDDE plus SED-sequence 
blazar model here, there are more hard-spectrum sources 
with lower gamma-ray flux. This trend already can be 
seen in the plotted blazar points in Fig. [T] 



C. Blazar Measurements 

Through one year of running, the Fermi-LAT has de- 
tected a total of 296 FSRQs, 300 BL Lacs, and 72 blazars 
of unknown type. The observed FSRQs have an average 
spectrum with photon index 2.48 and BL Lacs have av- 
erage photon index of 2.07 [38|. This power-law index is 
similar to the DGRB power-law index of 2.41, which sug- 
gests that unresolved blazars could be the primary source 
of the DGRB. Additionally, the stacked spectra of known 
blazars detected by the Fermi-LAT are responsible for 
15% of their total observed gamma-ray emission observed 
by the Fermi-LAT. The number of blazars observed above 
a given flux tends to follow a broken power law, with a 
break at F(> 100 MeV) = 6 x 10"^ photons cm'^ s"^ 
This break seems to be independent of detector sensitiv- 
ity, because the sensitivity dies off much more quickly as 
a function of flux than the blazar number count (FBIO). 

In the Fermi-LAT measurements, FSRQs and BL Lacs 
have similar variability properties, so the assumption 
that they are of one class appears valid. For BL Lacs, 
the LAT has detected significantly more hard-spectrum 



sources than soft-spectrum sources, which is consistent 
with the known selection bias in the measurement. FS- 
RQs peak at a redshift of unity, indicating that the sam- 
ple is approaching completeness. In contrast, BL Lacs 
peak at low redshift, indicating that the sample is not 
yet complete. FSRQs tend to be more luminous than 
BL Lacs: FSRQs have radio luminosities that peak at 
-^^rad ~ 10^^'^ erg/s whereas BL Lacs have lower radio 
luminosities peaking at Liad ~ 10''^ erg/s [35]. This 
would indicate that there is a fairly large contribution 
of low-luminosity, soft-spectrum BL Lacs that has yet to 
be resolved. 

The differences in spectra between FSRQs and BL Lacs 
are significant. The average gamma-ray photon index is 
roughly 0.5 larger for FSRQs than for BL Lacs. Even 
among BL Lacs themselves, high-synchrotron-peak BL 
Lacs have a photon index of 2.28 while low-synchrotron- 
peak BL Lacs have a photon index of 1.96. FSRQs give 
off their peak synchrotron radiation at around 10^'^ Hz 
whereas for BL Lacs, the distribution is much broader, 
stretching from lO^^ Hz to lO" Hz [38]. FSRQs have 
their inverse Compton (IC) peaks at energies less than 
100 MeV, so power-law fits work fairly well to match 
their LAT-measured spectra. For BL Lacs, the peak IC 
emission tends to lie in the LAT's energy range, with 
low-synchrotron-peak BL Lacs peaking closer to 100 MeV 
and high-synchrotron-peak BL Lacs peaking closer to 100 
GeV. Because of these peaks, these spectra do not match 
a power-law, though a broken power-law can approxi- 
mately fit them [28) . 

To truly model the blazar SED, a multiwavelength 
analysis is needed [27]. The Fermi-LAT Collaboration 
did a multiwavelength study of the spectra of blazars, 
combining the results of several radio, x-ray, optical, and 
gamma-ray blazar studies [3^]. This study found strong 
correlation between the x-ray and gamma-ray spectral 
slopes, indicating that blazar spectra fit a two-peaked, 
synchrotron plus IC scenario well. They found that BL 
Lacs have larger synchrotron peaks than FSRQs, which 
explains why BL Lacs have harder gamma-ray indices. 
This study plotted the SED for several blazars, all of 
which have a strong double-peaked shape when luminos- 
ity is plotted versus frequency on a log- log plot. This is 
consistent with previous analyses of the blazar SED [27] . 



III. DETERMINATION OF BLAZAR FLUX 
AND SPECTRUM 

A. Spectral Energy Distribution 

The model of blazar emission we use consists of two 
parts: a GLF to give the density of blazars per unit 
luminosity and an SED to determine the luminosity of 
blazars as a function of energy. These are denoted 
by p^{L^,z) and vL^{x;P), where z is redshift of the 
blazar, is the gamma-ray luminosity (defined as vL^, 
at hiy — 100 MeV), x = log]^Q(i^/Hz) for blazar rest- 
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frame frequency and P is the bolometric luminosity. 
Because our SED separates blazars according to radio lu- 
minosity, the bolometric luminosity is used to determine 
which SED curve matches a given blazar. For a given 
SED curve, the bolometric luminosity can be calculated 
as / Lj^dv. This can then be used to find the gamma-ray 
luminosity. 

Ref. ^27) analyzed the relationship between frequency 
and luminosity for blazars. To get these relationships, 
blazars were binned by radio luminosity. This analy- 
sis showed that blazar gamma-ray index is correlated 
with blazar luminosity. This correlation is consistent 
with the experimental results that FSRQs have high lu- 
minosities and large gamma-ray spectral indices while 
BL Lacs have lower luminosities and smaller spectral in- 
dices [551 [Ml HO] ■ A proper calculation using blazar spec- 
tra should account for this relationship between index 
and luminosity, and not simply use a power law in energy 
for the blazar spectrum. Note that this was not done in 
Ref. [1] , which claimed that blazars cannot constitute the 
full DGRB. 

For the frequency dependence of the blazar luminos- 
ity, we use the SED sequence of Inoue and Totani [8]. 
In this model, blazars SEDs are fit over frequencies from 
radio to gamma ray, as in Ref. 27\. Each SED is com- 
prised of two components, a synchrotron component at 
lower energies and an IC component at higher energies. 
These are each parameterized by a parabolic peak with 
a lower-energy linear tail. The details of the model are 
determined by fitting to the data in Ref. (27i , which give 
vLi, as a function of rest-frame frequency v for five lu- 
minosity bins. This provides the gamma-ray luminosity 
{vL^ at hv = 100 MeV), the specific luminosity L^{y)^ 
and the bolometric luminosity / L^dv for a blazar with 
known radio band luminosity (z^L^ at 5 GHz). The full 
model can be found in Appendix \K\ 

As a check on the versatility of the SED model, we 
explicitly compared the model to several blazar spectra 
measured by the Fermi-LAT Collaboration [28l|39]. The 
model fit the data in the Fermi-LAT energy range well. 
It also matched the data qualitatively: the model spec- 
tra had increasing, decreasing, or flat spectral shapes in 
agreement with the Fermi-LAT-measured spectra. Such 
agreement indicates that this SED fit approximates the 
full blazar SED well. 



Reference [5] argues that the gamma-ray luminosity 
can be related to the x-ray AGN disk luminosity Lx 
through the bolometric luminosity by P = lO'Lxi where 
q is a scaling parameter. This is because the bolometric 
luminosity from a blazar jet is proportional to the mass 
accretion rate m. For blazars with low accretion rate, the 
conversion of power into luminosity is inefficient, with 
Lx TV? . For blazars with high accretion rate close to 
the Eddington limit, the conversion is efficient and the 
disk luminosity goes as Lx oc m [42ti44] . Because black 
hole growth takes place mostly near the Eddington limit, 
it is reasonable to assume that P oc m oc Lx 05]. Note, 
Lx is the x-ray luminosity from the accretion disk of the 
blazar, not to be confused with the x-ray luminosity of 
the beam. 

The comoving number density per unit L-y of gamma- 
ray blazars is 

Pi{L^,z) = k'^^Px{Lx,z), (3.1) 

where px is the comoving number density of AGN per 
unit Lx, z is the redshift to the source, and k is the 
fraction of AGN observed as blazars. The quantity is 
referred to as the GLF. A parameterization of the x-ray 
luminosity function px is found in AppendixjB] The GLF 
has three free parameters: q determines the ratio of bolo- 
metric jet luminosity to accretion-disk x-ray luminosity, 
7i is the faint-end index that determines how the GLF 
behaves for low luminosities, and the blazar fraction is k. 

These GLF models are based on LDDE of AGN, as op- 
posed to PLE models. In PLE models, AGN luminosity 
changes with redshift, but the comoving density of AGN 
remains constant. This has been a popular method of 
determining blazar parameters [IHlEQl EHESJ • LDDE 
models have a peak evolution redshift which depends on 
luminosity, so AGN of different luminosities will have 
slightly different evolutions [29l [41] . This gives a better 
fit to the AGN data and should describe blazar evolution 
more fully than PLE models [23j. The exact relation- 
ship between x-ray AGN and gamma-rays blazars is not 
yet known. We are using the simple ansatz that they 



are related as shown in Eq. (3.1), as proposed by In- 
oue and Totani [5]. To the best of our knowledge, this 
model satisfies all current observations and constraints, 
and therefore is a viable possibility. 



B. Gamma-ray Luminosity Function 

For the distribution of gamma-ray blazars, we fol- 
low the hard x-ray AGN distributions parameterized by 
Ueda et al |29j . Similar work was done for soft x-rays 
by Hasinger et al [H]. However, the hard x-ray pa- 
rameterization gives a more conservative prediction of 
blazar detection by Fermi-LAT, so we use that here. 
For rest-frame (emission frame) energy of Egam.ros = 
100 MeV, the gamma-ray luminosity is given by L^ = 

(^gani,res/^)-^i/ (^gam,rcs/^; -^)' 



C. Calculation of Blazar Number and Flux 

For a given blazar, the gamma-ray flux observed on 
Earth is 

F,{z,P) = -i-ti- r (3-2) 

where is the luminosity distance, P is the bolometric 
luminosity, and i?min,obs — 100 MeV is the minimum 
observable photon energy on Earth by the Fermi-LAT. 
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With the GLF and SED, the number count of blazars 
detected above a sensitivity is 



iV(> F^) 47r 



dz- 



dV 



(3.3) 

where ilj™ is the luminosity below which a blazar at 
redshift z is no longer detectable for the sensitivity F^. 
We set the parameter Zmax = 5, but this does not affect 
the calculation significantly, since the peak distribution 
is at redshift of order unity. 

The diffuse flux coming from unresolved blazars is 
given by 



dN 



1 



dEjodAdtd^ An Jq 



dz 



(3.4) 



Here, E^ is the emitted photon energy [and E^^ — 
E.^/{1 + z) is the observed photon energy at Earth], A is 
area on Earth, t is time on Earth, and $7 is solid angle in 
the sky. Here L^/{E^) is the number of photons emitted 
per rest-frame frequency per rest-frame time per blazar 
(h is Planck's constant). The quantity dL-^p^ is the num- 
ber of blazars per comoving volume. The integral dx is 
the line-of-sight integral over the comoving distance. Be- 
cause for 7i > 1 the integral diverges at zero luminosity, 
i7,min is a lowcr bound on the luminosity integral. We 
choose L-y,min = 10^^ erg s~^ which is an order of magni- 
tude lower than any Fermi-LAT observed blazar [551 EH] ■ 
That is, we impose a step-function cutoff of blazar GLF. 
The final result is not strongly dependent on the value 
of this cutoff, with a 2-order-of-magnitude difference in 
F7,min modifying our best- fit parameters by ^25%. 

The exp(— r) factor in the diffuse flux calculation ac- 
counts for absorption of the photons on intergalactic 
background radiation before reaching Earth. We use the 
absorption factor from Gilmore et al. [46]. This absorp- 
tion factor was determined through the use of galaxy 
formation models to find the contribution of starlight to 



the absorption, as well as a contribution from quasars 
which is calculated based on empirical data. This model 
predicts lower values of the opacity t than previous es- 
timates, which leads to less expected absorption. This 
is consistent with the Fermi-LAT observing several high- 
energy photons coming from fairly high redshifts, and 
this opacity is consistent with the findings of Ref . [47] . 



D. DGRB Spectrum Calculation 

In addition to the blazar contribution to the DGRB 
flux, we also include a nonblazar AGN component to our 
DGRB spectrum calculation. Ref. [lOJ has shown that 
nonblazar AGN can account for the background radiation 
down to keV energies. The combination of blazars with 
nonblazar AGN gives a unified model that can explain 
the diffuse high-energy x-ray to gamma-ray background 
over 8 orders of magnitude in energy. 

The AGN model we use is the model of Ref. lOj. This 
model assumes the usual thermal electrons from AGN 
coronae, but it includes a high-energy nonthermal com- 
ponent as well. These electrons Comptonize, which pro- 
duces the known x-ray spectra of AGN. This high-energy 
component is analogous to the emission from solar coro- 
nae in solar flares. Such electrons are assumed to have a 
power-law injection spectrum dN/dE cx E~^ . By adding 
this nonthermal electron source to the usual thermal one, 
it is found that the model matches the diffuse background 
spectrum well from energies from keV to tens of MeV. 

Specifically, we choose the F = 3.5 nonblazar AGN 
model of Ref. [10] , which we increase in amplitude by a 
factor of 2 in order to match 50% of the amplitude of the 
lowest-energy point in the Fermi-LAT DGRB spectrum, 
with a broken power law matching the measurements of 
the diffuse background by the COMPTEL Collaboration 
[48) . The power-law slope of the nonblazar AGN spec- 
trum is fixed by modeling of the hard x-ray luminosity 
function from x-ray AGN surveys [TUl US] , and the ampli- 
tude is fixed to match the lowest point in the Fermi-LAT 
DGRB spectrum. This amplitude is fixed throughout our 
fitting. In order to refiect the uncertainty of the ampli- 
tude of the flux in the lowest-energy bin, we allow for 
it to have an amplitude uncertainty of 10%, which we 
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vary and show in Fig. [Sj Another low-energy emission 
source, such as miUisecond pulsars or star-forming galax- 
ies, may be responsible for the lowest-energy portion of 
the DGRB, but our analysis is not strongly dependent 
on the spectral shape taken by the low-energy emission 
source. For example, the gamma-ray spectrum from star- 
forming galaxies in Ref. [11] has a similar shape and po- 
tential amplitude as the nonblazar AGN component. 

In our blazar model, there are three free parameters, 
in addition to those fixed in the nonblazar AGN model, 
as described in Sec. |IIIB| q, 71, and n. All other parame- 
ters in the blazar model are fixed to values based on data 
from other observations such as the SED sequence. It is 
the purpose of this paper to determine how well unre- 
solved blazars can reproduce the DGRB. Therefore, we 
simultaneously fit to the blazar source-count distribution 
dN/dF from Ref. [4| and the DGRB spectrum from FSIO. 
This simultaneous fit allows some freedom in the blazar 
spectrum while still conforming to known blazar num- 
ber distributions. We can use the results of such a fit to 
constrain models of the DGRB from unresolved blazars 
and predict a consistent model of the 5- year Fermi-LAT 
measurements of the DGRB. 

Fitting the model to the blazar dN/dF and the DGRB 
spectrum, we found that a simultaneous fit was quite 
reasonable. We set the lowest blazar luminosity as 
-^7,mm = 10"*^ erg s~^, as discussed above. The best- 
fit values we get are q = A.WtoH, 7i = l-Sllo.Jg, 
and logio(«;/10-'^) = 0.38l[J;7^ (95% CL). The best-fit 
68% and 95% CL regions for q and 71 are shown in 
Fig. [2) These are consistent with previous work [5], 
though more constrained because we are also fitting the 
source-count distribution function dN/dF. The model 
reproduces the DGRB and blazar dN/dF, with a reduced 
X^/DOF = 0.63. The value of q indicates that the bolo- 
metric luminosity of a blazar jet is roughly 15 thousand 
times more luminous than the x-ray from the accretion 
disk. Here, 71 > 1.0 so low- luminosity blazars have sig- 
nificant contributions to the total blazar fiux. There- 
fore, a ten or more order-of-magnitude lower value of 
^7,min would modify the calculation considerably, though 
no blazars have been detected below our i^^min thresh- 
old, and therefore it seems unlikely that there is a large 
population of very-low-luminosity blazars. The fraction 
K ~ 2.4 X 10~^ implies that there is roughly one blazar 
for every 420 thousand nonblazar AGN. Our fit to the 
DGRB spectrum is shown in Fig. [3] and the fit to dN/dF 
is in Fig. |4] 

Our value for the AGN XLF and blazar GLF ratio k, 
3.4 X 10-6 to 5 X lO"'^ (at 95% CL), is similar to and 
slightly larger than the central value derived by Inoue & 
Totani [8], 1.7x IQ-^. This implies that only a small frac- 
tion of x-ray loud AGN is visible as gamma-ray blazars. 
The intrinsic jet opening angle of a blazar has been found 
to be ~ 1 deg (subtending an area of ~ 2 x 10"^ stera- 
dian) [5T]. Following from this is that only ^ 2 x 10^^ 
of the AGN jets are potentially visible as blazars. Our 
model then requires that only < 20% of AGN jets are 



gamma-ray blazars. This is not inconsistent with jet 
models [52], though if this fraction drops considerably 
(i.e., K is required to be much smaller), then it would call 
into question the blazar model analyzed here. 

Note that using the dN/dF estimated from a power- 
law blazar spectrum model is not perfect, due to the 
fact that the detection efficiency estimate depends on the 
spectral model [1]. However, Ref. 4j tested the dN/dF 
dependence on the sensitivity estimate with a non-power- 
law fit to the blazar spectra and found it did not sig- 
nificantly change the measurement of dN/dF. We also 
verified this sensitivity dependence with a test fitting by 
increasing the errors on the measured dN/dF at low flux, 
and we found that our model did not prefer a different 
amplitude or shape to the source counts at the low flux 
where the efficiency for blazar detection is low. 

Refs. [3 [9] used a combined GLF plus SED model to 
predict the Fermi-LAT's ability to observe blazars and 
their spectra, using the results of the EGRET Collabo- 
ration. The paper fit its GLF parameters using the red- 
shift and gamma-ray luminosity distributions of EGRET 
blazars. This led to a prediction that 600 to 1200 blazars 
should be resolved in 5 years of Fermi-LAT data, which 
would yield 98% to 100% of the total blazar flux. How- 
ever, the cumulative number of blazars predicted by that 
paper is in disagreement with the observations of the 
Fermi-LAT W . The cumulative number count by Ref. [5] 
is predicted to have a break at 10^^ photons cm~^ s~^ 
whereas the break seen by the Fermi-LAT Collaboration 
is at 5 X 10~* photons cm~^ s~^. Also, the surface den- 
sity of sources predicted in that paper is too small to 
match the measured value. 

Importantly, Refs. fit their model to the EGRET 
catalog blazar spectra SED, not that from Fermi-LAT. 
The EGRET telescope had strong cuts which limited 
high-energy photon observations, which lead to EGRET 
only observing a few BL Lacs [4]. Also, the redshift and 
luminosity distributions are strongly dependent on de- 
tector sensitivity, because BL Lacs have lower luminosity 
and therefore are observed at lower redshifts. This means 
that the current data for the overall blazar redshift dis- 
tribution, in particular, is more strongly biased toward 
lower redshifts than the complete distribution. Ref. [S3] 
posited that one significant source for the difference be- 
tween this calculation and the Fermi-LAT results comes 
from needing to correctly account for Fermi-LAT sensi- 
tivities. By fitting to dN/dF, which is not as heavily 
dependent on detector sensitivity, we can get a more ro- 
bust prediction that should not change significantly for 
different sensitivities. Refs. [8j 9J argued that a model of 
this type should roughly match the DGRB spectrum. In 
Ref. f8], the model parameters were fit to the EGRET 
DGRB spectrum, and, as discussed above, the model pa- 
rameters are roughly consistent with our results. In our 
analysis here, we use the DGRB spectrum and fiux source 
counts, as measured by the Fermi-LAT, as a constraint 
in order to determine how well this class of models fits 
the DGRB and blazar population. For those models that 
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FIG. 3. Shown are the best-fit model for the current DGRB spectrum (sofid black line) and our upper/lower 95% CL forecast 
for the Fermi- LAT 5-year sensitivity (magenta-star/green-circle points). The low-energy-dominating solid red line is the AGN 
flux from Ref. [10) . The high-energy-dominating blue lines are the blazar contribution to the DGRB for the current (solid), and 
predictions for the most-optimistic (dashed) and least-optimistic (dotted) 95% CL 5-year Fermi-LAT resolved fractions. The 
grey lines are the combined 95% CL AGN plus blazar predicted flux for the corresponding blazar contribution. The DGRB 
data (triangles) are from FSIO and the COMPTEL data (diamonds) are from Ref. [50] . 



fit the spectrum, we can determine the predicted values 
for the DGRB flux at the Fermi-LAT's 5-year sensitivi- 
ties and determine the theoretical uncertainty on these 
predictions. 

In another analysis of the contribution of blazars to the 
DGRB, the Fermi-LAT Collaboration used the currently 
measured differential number distributions of blazars 
{dN/dF) and blazar gamma-ray index (F) distributions 
to estimate the contribution of unresolved blazars to the 
DGRB [1]. In that analysis, it was found that less than 
20% of the DGRB can be accounted for by blazar emis- 
sion. However, in that calculation, the assumption was 
made that the distribution of indices F is independent 
of sensitivity. Because less-luminous BL Lacs have sig- 
nificantly different indices than more luminous FSRQs, 
the overall distribution of indices should change as bet- 
ter sensitivity allows a greater fraction of BL Lacs to be 
detected. 

Additionally, it was shown in Refs. [32] that a ba- 
sic power-law model does not fit the individual blazar 
spectra well, especially for the low-luminosity BL Lacs. 
A GLF plus SED model should overcome these issues. 
The GLF accounts for differing redshifts of blazars, so 
the relationship between flux sensitivity and luminosity 



detectability is well-defined. The SED accounts for the 
distribution of luminosities with energy, so a calculation 
around the IC peaks for BL Lacs should more realistically 
reproduce the contribution to the DGRB from blazars 
than a simple distribution of photon indices. This is es- 
pecially important to incorporate when determining the 
contribution of unresolved low-luminosity blazars to the 
DGRB, since they have much harder spectra than high- 
luminosity blazars. 



IV. 5- YEAR PREDICTIONS FOR BLAZARS 
AND THE DGRB 

We adopt the 5-year predictions for a sensitiv- 
ity to point-sources by Fermi-LAT of ^5 ~ 2 x 
10~^ photons cm~^ s~^ above 100 MeV. This value is 
consistent with the Fermi-LAT Collaboration's estimate 
of the LAT sensitivity to point sources with gamma-ray 
index of ~2 [53] P] As discussed earlier, the majority of 



^ http : //f ermi . gsf c . nasa . gov/science/433-SRD-0001_CH-04 . pdf 
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FIG. 4. Shown is the best-fit model for the source-count dis- 
tribution function dN/dF (solid line). The data are from 
Ref. [H 




■^100 s ^ cm ^) 



FIG. 5. Shown is the flux per logarithmic sensitivity for 
our best-fit model. The dashed line is the flux coming from 
blazars and the solid line is the flux coming from nonblazar 
AGN. The vertical solid lines with arrows mark the sensitivity 
to all spectral-index sources at the Fermi-LAT 1-year and the 
projected 5-year sensitivity of Fermi-LAT. The gray boxed 
region indicates the range of sensitivity at 1-year to sources 
with the spectral indices of the bulk of the blazar population, 
as in Fig. [l] 



low-flux blazars are expected to be BL-Lacs, which pre- 
dominantly have radio luminosity less than 10^'^ erg/s 
[55] , Such low- luminosity blazars have gamma-ray in- 
dices of ~2 or less, according to the blazar SED. There- 
fore, we find the use of ^5 = 2x 10~^ photons cm~^ s~^ as 
the Fermi-LAT 5-year sensitivity to blazars of all gamma- 
ray indices to be a reasonable estimate. 




40 42 44 
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FIG. 6. Shown is the radio luminosity distribution of blazars. 
The solid line is our prediction for the distribution after 5 
years of Fermi-LAT running. The dotted line is the current 
Fermi-LAT distribution for blazars [38 . Each distribution is 
independently normalized to unity. 




FIG. 7. Shown are the distribution in redshift of blazars. 
The solid line is our prediction for the distribution after 5 
years of LAT running. The dotted line is the current Fermi- 
LAT-measured distribution for FSRQs and the dot-dashed 
line is the current Fermi-LAT-measured distribution for BL 
Lacs [38]. Each distribution is independently normalized to 
unity. 



To determine the total number of blazars detectable 



by the Fermi-LAT, we need to take Eq. (3.3) down to a 
sensitivity of 6*5. Similarly, we can determine the total 
number of blazars in the sky by letting the sensitivity go 
to zero fiux. With 95% CL, we predict that there are 
5.4^j^ y X 10^ total blazars in the observable universe. Of 
these, 2415^420 should be detectable by the Fermi-LAT 
after 5 years of running. The amount of flux coming from 
blazars per logarithmic sensitivity is shown in Fig. [Sj Our 
prediction is that 94.7^2''i9o of blazar flux is expected to 
be resolved by the Fermi-LAT after 5 years, mostly at 
lower energies. In contrast, the flux for nonblazar AGN 
should not be appreciably resolved for another 4 orders 
of magnitude in sensitivity. 

In addition to the number counts of blazars, we can 
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also predict the distributions of blazars in luminosity 
and redshift. To get these distributions, we differenti- 



ate Eq. (3.3). The distribution of blazars in radio lumi- 



nosity, shown in Fig. [6] shifts toward lower luminosities 
at better sensitivities. This is due to the FSRQ popu- 
lation being mostly resolved, whereas the new resolved 
sources at better sensitivities are mostly low-luminosity 
BL Lacs. The redshift distribution of blazars, Fig. [7j 
should shift toward higher redshifts as sensitivity im- 
proves. Because the FSRQ sample is mostly complete, 
it would be expected that the redshift distribution of BL 
Lacs, and blazars in general, should be roughly similar to 
the current redshift distribution of FSRQs. Our predic- 
tion of the redshift distribution of blazars after 5 years of 
Fermi-LAT running matches well with the current FSRQ 
distribution, which provides a verification of our theory 
and fit parameters. Note that the FSRQ sample is not 
totally complete, and the objects to be resolved at z > 2 
would be FSRQs. As can be seen in Fig. 9 of Ref. [55] , 
the distribution of FSRQs reaches the current flux limit, 
so there remains a population of high-luminosity, soft 
spectral index, high redsift FSRQs to be resolved. 

In our model, we fit the total blazar plus AGN flux to 
the DGRB spectrum for the spectrally-dependent sen- 
sitivity as described above. The model fit worked ex- 
ceptionally well, indicating that a combination of blazar 
flux with the flux of nonblazar AGN makes up all the 
DGRB over a wide range in energies. With this fit, we 
then calculated what the combined flux should be after 
5 years of Fermi-LAT observations, giving the sensitiv- 
ity of 2 X 10^^ photons cm^^ s^^. The upper and lower 
bounds of the 95% CL region of this calculation are given 
by the upper and lower forecast points in Fig. [3] We 
have included a 10% uncertainty on the nonblazar AGN 
flux in this error estimate to account for the error in the 
lowest-energy bins' constraint on the AGN model. At 
100 GeV, we expect the DGRB to decrease by a factor 
of 1.6 to 2.6 at the 95% CL upper and lower flux limits, 
whereas at 100 MeV the DGRB only decreases by a fac- 
tor of 1.3 to 1.9. The difference in DGRB improvement 
is due to a greater fraction of the DGRB being due to 
blazars at high energies, while the nonblazar AGN flux 
dominates at low energies. Importantly, the resolution of 
sources can do better than the square root of exposure 
time due to the increased prevalence of easily-detected 
hard sources beyond, but near, the current point-source 
flux-limit sensitivity. 



V. CONCLUSIONS 

We have shown that the DGRB can be composed en- 
tirely by gamma-rays produced in blazars and nonblazar 
AGN. The LDDE plus SED-sequence is a physical model 
for the spectral evolution of a cosmologically-evolving 
blazar population contributing to the DGRB based on 
the unified AGN model for blazars. This model success- 
fully accounts for the full DGRB spectrum as well as 



the full blazar source-count distribution function, which, 
unlike other approaches, are not used as components of 
the model. Independent of the nonblazar AGN compo- 
nent, the blazar model produces nearly the entire DGRB 
at its highest measured energies. The small value of 
K ~ 2.4 X 10^^, the x-ray AGN fraction seen as blazars, 
constrains this model to require a small fraction, <20%, 
to be both properly oriented and sufficiently energetic in 
order to be gamma-ray emitters. 

We found constraints on this model from the spec- 
trum of the DGRB and source-count distribution func- 
tion dN/dF of blazars as observed by Fermi-LAT. Our 
results are consistent with previous work by Inoue & 
Totani [8 which employed EGRET spectral data to fore- 
cast the Fermi-LAT DGRB. We forecast that 94.7^2:?% 
of the fiux from blazars will be resolved into point sources 
by Fermi-LAT with 5 years of observation, with a corre- 
sponding reduction of the flux in the DGRB by a factor 
of ^2 to 3 (95% CL) from the automatic removal of these 
sources in the measurement of the DGRB. This has sig- 
nificant consequences for the sensitivity of the DGRB 
measurement to dark matter annihilation, which we ex- 
plore in a companion paper |55| . 

We predict that 2415^42° blazars should be resolved, 
of 5Atl[Y X 10* total blazars in the universe (95% CL). 
Recent results of anisotropy in the DGRB also indicate 
the likely presence of an unresolved point-source popu- 
lation |56] . Using tests with enhanced point-source sen- 
sitivity, we find that future gamma-ray experiments at 
Fermi-LAT energies will resolve the blazar contribution 
to the DGRB such that the flux in the DGRB decreases 
as the square root of the point-source sensitivity. 

The LDDE plus SED-sequence model is more com- 
plex than the over-simplistic source-count method with 
a fixed spectral-index distribution adopted by the Fermi- 
LAT Collaboration in FBIO, yet it has fewer free param- 
eters for the blazar population than the more simplified 
model (three versus four free for the blazar model, plus 
those fixed in the nonblazar AGN model in this work). 
Most importantly, the Fermi-LAT analysis of FBIO fixes 
the spectral index of the blazar population, and, cru- 
cially, does not include the hardening of the spectra of 
the unresolved low-luminosity blazar population. The 
hardening of spectra with lower luminisity has been seen 
by both EGRET [27] and Fermi-LAT (Fig.[l]). The fixed 
spectrum forces the FBIO conclusion that only ~16% of 
the GeV isotropic diffuse background could arise from 
blazars, and is also the case in other work using fixed 
blazar spectra _30J. Other recent work with different 
blazar population models, including spectral shape vari- 
ation [ST], possible point-source confusion [5S], and BL 
Lac dominance of the unresolved portion I59j also find 
that a substantial portion of the DGRB could arise from 
the blazar population. 

Overall, the SED-sequence model of blazars and AGN 
as the source of the DGRB is remarkably consistent with 
the measured DGRB spectrum and blazar source-count 
distribution. The SED-sequence will continue to be im- 
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proved with upcoming Fermi-LAT blazar data [BD] . Fur- 
ther analyses of the type presented here, incorporating 
potential enhancements to the SED-sequence model, the 
XLF of AGN, and general studies of observed blazar spec- 
tral properties, will further enlighten the understanding 
of the extragalactic gamma-ray sky. 
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Appendix A: Blazar SED Sequence 

The full SED fit is given as a function of radio lumi- 
nosity ipR and the logarithm of rest-frame frequency x. 
We follow Ref. in the formulation of the SED. The 
radio luminosity is used to distinguish between SEDs for 
blazars of different bolometric luminosity. This separa- 
tion of SED by total luminosity should account for the 
difference in spectral index seen by the Fermi-LAT be- 
tween higher-luminosity FSRQs and lower-luminosity BL 
Lacs [21131110]. 

x = logioWHz), (Al) 



ijj{x; tpR) = lo, 



erg s~ 
= i^{x = 9.698). 



(A2) 
(A3) 



The full model is the sum of a synchrotron ['(/'^(a;)] and 
inverse Compton [^/;c(a;)] component. 



V'(a;) ^logioilC^^^^) -K10'^<=(^)] 



(A4) 



Each component is parameterized as the sum of 
a lower-frequency linear part and a higher-frequency 
parabolic part. Here, Xtr,s and Xtr,c si's the frequencies 
where the linear part transitions to the parabolic part for 
the synchrotron and IC component. The linear parts are 
written as 

i)sl{x) = (1 - as){x - Xr) ^^I^r {x < Xtr,s), (A5) 

ipci{x) = (1 - ac)ix - xx) +ipx {x < xtr,c), (A6) 

where as — 0.2 and etc = 0.6 are the L^, oc indices 
in the radio and hard x-ray bands, respectively. The 
characteristic radio and hard x-ray frequencies are xr = 
9.698 and xx = 17.383. The radio luminosity 4'R is 



an input parameter to the theory and the hard x-ray 
luminosity is fitted to the data as 

r (^fl- 43) +43.17 Vii<43 

ipx ^ { 1.40(V'fi- 43) 43.17 43 < < 46.68 
[ 1.40(46.68 - 43) + 43.17 ^pR > 46.68 . 

(A7) 

The parameter V'x is kept constant for ipR > 46.68 be- 
cause the continuity of the IC component cannot be sat- 
isfied above this value. However, this hard x-ray lumi- 
nosity corresponds to a gamma-ray luminosity well above 
the maximum detected gamma-ray luminosity, so it does 
not affect the calculation of the DGRB. 

The parabolic parts of the components are parameter- 
ized as 

1ps2ix) = il^s.p - [{x - Xs)/a]'^ (x > Xfr^s), (A8) 
1pc2{x) = 1pc,p - [{X - Xc)l<jf {X > Xtr,c), (A9) 

where Xs and the synchrotron and IC peak fre- 

quencies, tps.p and ipc,p are the synchrotron and IC peak 
luminosities, and a is the width of the parabolas. 

By requiring continuity of the synchrotron component 
from the linear-to-parabolic parts, we have 



{l~as){xtr,s-XR)+'lJjR + 



(AlO) 



Similarly, the continuity of the IC component gives 

-*.c = '^'f^ , (All) 

C = a^{l-a,)-2xc. (A12) 

ri = xl + a^lijjx - xxil- ttc) - i'c.p] ■ (A13) 
By inspection 

xtr^s = 10.699, (A14) 

Xc^Xs + 8.699 . (A15) 
Fitting to data, the rest of the parameters are given by 
-0.88(V'fl - 43) + 14.47 ipR < 43 



-0.40(V'i? - 43) + 14.47 ipR > 43 
0.0891X, + 1.78 Vfl < 43 

[2{xs-xtr^s)/{l-as)f^ ijR>43 



(A16) 
(A17) 



-Ips^p 



V'fl, < 43 



1.77(V'fl-43)°-™+45.3 > 43 ^^^^^ 



These parameters have been chosen such that the lumi- 
nosity changes continuously with tJjr over all luminosities 
and to make the synchrotron linear-to-parabolic transi- 
tion smooth for large iJjr. 



Appendix B: X-ray Luminosity Function 

The x-ray luminosity function px is the comoving num- 
ber density of AGN per unit x-ray AGN disk luminosity 
Lx- The model of Refs. [HUT] give the distribution as 



px{Lx,z) = pxiLx,0)f{Lx,z). 



(Bl) 
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The present distribution is given by 



Px{Lx,Q) 



Ax 
Lxln{lO) 



Lx 



71 



r \ 721 -1 



The density evolution is given by 

r (i + 2)pi 

f{Lx,z) = 



(B2) 



l+z 



The peak evolution happens at z^, given by 

(T \ _ j ^x > La 

^-(^^)-\ z*{Lx/Lar Lx<La. 

The evolution indices pi and p2 are 



pi +/3i[logio(ix)-44.0] 
P2=p;+/32[logio(ix)-44.0]. 

The parameters for the models are given in Table I. If 
7i > 1, then the integrated background flux diverges, so 



we set the minimum gamma-ray luminosity to L-y.min = 
10^^ erg/s. This is an order- of- magnitude lower than 
any Fermi-LAT observed blazar, and the results are not 
sensitive to this value being lowered slightly [251 132] • 

Table 1 



Parameters for the AGN XLF 



ciLx) 


Parameter 


Ueda et al. 2003 


Hasingcr ct al. 2005 


ciLx). 
(B3) 


Ax (Mpc-3) 
logioL\ 

72 


5.04 X 10"*^ 

A'X 04+0 21 

2.23 ±0.13 


2.62 X 10-^ 
43.94 ±0.11 
2.57 ±0.16 




4 


1.9, fixed 


1.96 ±0.15 


(B4) 


logioLa 


44.6, fixed 


44.67, fixed 


a 


0.335 ±0.07 


0.21 ±0.04 




Pi 


4.23 ± 0.39 


4.7 ±0.3 




Pi 


-1.5, fixed 


-1.5 ±0.7 


(B5) 
(B6) 


Pi 


0.0, fixed 


0.7 ±0.3 


/32 


0.0, fixed 


0.6 ±0.8 



Note: Luminosities are in erg/s 
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